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In this work, the flame characteristics of torrefied biomass were studied numerically under high-temper- 
ature air conditions to further understand the combustion performances of biomass. Three torrefied bio¬ 
masses were prepared with different torrefaction degrees after by releasing 10%, 20%, and 30% of volatile 
matter on a dry basis and characterized in laboratory with standard and high heating rate analyses. The 
effects of the torrefaction degree, oxygen concentration, transport air velocity, and particle size on the 
flame position, flame shape, and peak temperature are discussed based on both direct measurements 
in a laboratory-scale furnace and CFD simulations. The results primarily showed that the enhanced drag 
force on the biomass particles caused a late release of volatile matter and resulted in a delay in the igni¬ 
tion of the fuel-air mixture, and the maximum flame diameter was mainly affected by the volatile con¬ 
tent of the biomass materials. Furthermore, oxidizers with lower oxygen concentrations always resulted 
in a larger flame volume, a lower peak flame temperature and a lower NO emission. Finally, a longer 
flame was found when the transport air velocity was lower, and the flame front gradually moved to 
the furnace exit as the particle size increased. The results could be used as references for designing a 
new biomass combustion chamber or switching an existing coal-fired boiler to the combustion of 
biomass. 

© 2013 The Combustion Institute. Published by Elsevier Inc. All rights reserved. 


1. Introduction 

It has been widely recognized that biomass has the potential to 
fulfill the increasing energy demands and the strict regulations on 
C0 2 emissions m. Sweden has set a target to create the first oil- 
free economy in the world by 2020, and today, approximately 
21% of the total heat of consumption in Sweden is provided by bio¬ 
mass and biofuels [2]. Currently, the combustion of biomass is 
acknowledged to be a cheaper, lower risk, and more technically 
complete option compared to other biomass thermal conversion 
technologies. The combustion behavior of biomass is different from 
that of coal. Firstly, biomass yields a much higher fraction of vola¬ 
tiles through its devolatilization process. In addition, the low den¬ 
sity of biomass results in the higher oxidization rate of biomass 
char. Therefore, although a lower replacement ratio of biomass in 
co-firing boilers shows relatively few operational problems [3], 
the co-firing of greater than 20% of biomass (on an energy basis) 
in existing large-scale pulverized coal boilers without modifica¬ 
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tions usually results in unstable boiler efficiency and fluctuations 
of the boiler load. 

Torrefied biomass can be a suitable alternative fuel in existing 
large-scale pulverized coal boilers, because torrefaction changes 
the characteristics of biomass and moves it more similar to coal 
[4-9]. For example, torrefied biomass has higher energy content 
when compared to the parent material, the lower fibrous content 
makes it easy and energy-efficient to grind the torrefied materials, 
torrefied biomass powder have high flowability and good fluidiza¬ 
tion behavior and biomass is increasingly uniform after torrefac¬ 
tion. Accordingly, a torrefaction-based co-firing system in 
pulverized coal boilers has been proposed as a promising option 
with the aim of 100% of biomass fuel switching [10,11]. 

When co-firing torrefied biomass in an existing pulverized coal 
boiler, the biomass flame characteristics are important for select¬ 
ing the burner for biomass combustion from existing burners, 
and also may provide an operational reference for fuel switching. 
Several previous works have attempted to study the differences 
between the combustion of pulverized coal and biomass, or the ef¬ 
fects of co-firing on flame characteristics and stability. Lu et al. [3] 
suggested that the addition of biomass affected the flame charac¬ 
teristics, particularly the flame ignition points and brightness, 
due to the variation in the physical and chemical properties of 
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the biomass fuels. Ballester et al. [12] compared pulverized coal 
flames and pulverized sawdust flames in a semi-industrial furnace 
under the same operating conditions and concluded that the com¬ 
bustion zone is extended further for the biomass. Molcan et al. [13] 
found that the addition of biomass to coal could improve the com¬ 
bustion efficiency because of the lower CO concentrations and 
higher char burnout levels in co-firing. 

Computational fluid dynamics (CFD) has also been widely em¬ 
ployed to study biomass combustion by compiling a comprehen¬ 
sive biomass combustion model. Mando et al. [14] performed a 
CFD simulation of pulverized coal and straw combustion by using 
a commercial multi-fuel burner and concluded that straw combus¬ 
tion is associated with a significantly longer flame and smaller 
recirculation zones compared to coal. Yin et al. [15] numerically 
investigated the co-firing of a biomass suspension by studying 
the coal flame and the biomass flame in a swirl-stabilized burner 
and reported that a greater amount of volatiles released from the 
straw resulted in a greater volume in the reactor prior to mixing 
with the oxidizer. Gubba et al. studied the influence of particle 
shape and the internal thermal gradients of the biomass particles 
on the co-fired flames of pulverized coal and biomass [16]. Addi¬ 
tionally, other CFD simulations have been performed on large-scale 
furnaces to investigate the effects of the co-firing ratio on boiler 
performances [10,17-19]. However, there has not been a compre¬ 
hensive study on the effects of the torrefaction degree, fuel injec¬ 
tion velocity, oxygen concentration, and particle size on biomass 
flame characteristics with the purpose of applying torrefied bio¬ 
mass in a large-scale furnace either for pure biomass combustion 
or for co-firing. The CFD modeling of biomass combustion still 
faces significant challenges because of the diversity of the biomass 
fuels and the lack of knowledge regarding their combustion prop¬ 
erties [20]. In this work, the swelling behavior during devolatiliza¬ 
tion and the shape of torrefied biomasses have been 
experimentally studied and considered to accurately predict the 
particle trajectory and thus accurately predict the biomass flame 
and combustion properties. 

Moreover, volumetric combustion has been proposed and 
developed to realize a high replacement ratio of biomass in co-fir¬ 
ing boilers, which is expectedly accompanied by a reduction in pol¬ 
lutants [21]. The volumetric combustion of both gas and liquid 
fuels has been successfully realized for conditions of High Temper¬ 
ature Air Combustion (HTAC) [22-25]. Currently, several attempts 
have also been made to apply such technology in pulverized coal 
boilers [26,27]. Such a volumetric combustion could be character¬ 
ized as an extension of the HTAC technology to solid fuels. In this 
work, to further understand the combustion performance when 
switching a pulverized coal-fired boiler to combust biomass, the 
combustion and flame characteristics of torrefied biomass have 
been studied numerically in a laboratory-scale furnace at high 
temperature air conditions. Three torrefied biomasses and a pri¬ 
mary coal were studied and a total of 17 cases were simulated with 
CFD, in which coal was employed for comparison and validation. 
Finally, the effects of the torrefaction degree, oxygen concentra¬ 
tion, transport air velocity, particle size on the flame position, 
shape, and peak temperature are discussed. 


2. Methods 

2.2. Fuel parameters and kinetic studies 

Three torrefied biomasses (TB) were prepared in previous, and 
were obtained from the same parent biomass, palm kernel shells 
(PKS), with different torrefaction degrees after releasing 10%, 20% 
and 30% of the volatile matter on a dry basis. For simplification, 
these biomasses are abbreviated as TB1, TB2, and TB3, respectively. 


The proximate and ultimate analyses data of the torrefied biomass 
samples and coal are shown in Table 1. It can be observed that tor- 
refaction significantly reduces the volatile-to-fixed carbon ratio of 
the biomass and upgrades its energy density. 

The Van Krevelen diagram provides the relationship between 0/ 
C and H/C for different solid fuels, as shown in Fig. 1. Three torr¬ 
efied biomasses and a coal from this work are shown in this graph, 
while a bituminous coal, a sub-bituminous coal, a wood and the 
raw PKS are provided for reference. Figure 1 shows that the raw 
biomass has a lower energy density than coal; however, it interest¬ 
ingly shows that the energy density of the torrefied PKS becomes 
closer to that of bituminous coal as the torrefaction degree in¬ 
creases. Meanwhile, the breakdown of the fiber structure during 
torrefaction allowed the torrefied PKS samples to be grinded by 
the existing coal mill system. The coal particle size was assumed 
to follow the Rosin-Rammler distribution with an average diame¬ 
ter of 61.7 pm and spread parameter of 4.21, and the torrefied bio¬ 
mass particles were assumed to have the same sizes because they 
were grinded by the same mill system. 

To further understand the chemical properties of the torrefied 
biomass, thermogravimetric (TG) tests were performed in a TG bal¬ 
ance. A TA-Q500 analyzer was used in constant heating rate tests 
under nitrogen (to temperature of 800 °C) for studying the pyroly¬ 
sis step, followed by the oxidation in air for measuring the ash con¬ 
tent. Furthermore, the devolatilization tests were carried out in an 
isothermal plug flow reactor (IPFR). The test conditions were set to 
be similar to industrial pulverized fuel applications, such as high 
temperature and high heating rates. The reactor used can represent 
the conditions found in full-scale combustion processes with heat¬ 
ing rates on the order of 10 4o C/s, a maximum temperature of 
1400 °C, and particle residence time from 5 to 1500 ms. A flow of 
lOOg/h of biomass powder guaranteed stationary conditions for 
a time test of at least 30 min. The fuel conversion was evaluated 
by comparing the composition of the collected solid residue with 
that of the parent fuel by using the ash tracer method. 

After the devolatilization tests, a nominal temperature of 900 °C 
and residence time of 300 ms were used to produce a high heating- 
rate char and assure the complete devolatilization but prevent the 
graphitization of the carbonaceous matrix. The collected char res¬ 
idues of the torrefied biomass were tested in the ThermoGravimet- 
ric Balance to evaluate the intrinsic char oxidation kinetics. Small 
char samples (3 mg) have been tested to limit material and thermal 
resistances during the tests. Three trials were carried out for each 
char sample at a constant heating rate (20 °C/min) under different 
transport gas mixtures (5%, 12%, and 21% volume of oxygen, with 
the remainder nitrogen). In addition, the optical microscope 
images of a statistical number of particles were analyzed to obtain 
the morphological parameters of the shape factor and swelling 
coefficient, and all of the morphological studies referred to the pla¬ 
nar outlines of each particle. More detailed procedures of biomass 
characterization were reported by Biagini et al. [28]. 

2.2. Experimental methods 

The scheme for the pulverized fuel combustion facility is shown 
in Fig. 2. The secondary air was heated to 1150°C by a compact 
pre-heater and injected through the top of the chamber. The pul¬ 
verized fuel was pneumatically conveyed into the furnace by trans¬ 
port air from a 5 mm vertical water-cooling lance. On the bottom 
of chamber, there was outlet for the flue gas. The temperature pro¬ 
file along the combustion chamber and the temperature at the inlet 
and the outlet were recorded using thermocouples. The tips of the 
thermocouples were located in the middle of the inlet and outlet 
channels, and all of the thermocouples tips installed on the furnace 
chamber were 200 mm from the wall. Two isokinetic probes were 
located in the flue gas channel to collect ash and sample flue gas. 
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Table 1 

Analysis data of pulverized biomass powder. 


Items 

Units 

TB1 

TB2 

TB3 

Coal 

Torrefaction degree 

% 

10 

20 

30 

- 

Proximate analysis 

Moisture (as received) 

wt% 

5.8 

5.6 

5.1 

6.7 

Ash (dry) 

wt% 

6.56 

7.31 

8.5 

11.10 

Volatile (dry) 

wt% 

66.53 

61.87 

49.92 

19.41 

Fixed carbon (dry) 

wt% 

26.91 

30.82 

41.58 

62.79 

LHV (dry) 

MJ/kg 

20.5 

20.9 

25.0 

28.95 

Ultimate analysis (dry ashfree) 

C 

wt% 

54.30 

55.43 

60.67 

86.27 

H 

wt% 

5.07 

4.90 

4.33 

4.54 

0 

wt% 

33.52 

31.83 

25.90 

7.04 

N 

wt% 

0.50 

0.50 

0.57 

1.70 

S 

wt% 

0.05 

0.04 

0.03 

0.45 



Fig. 1. Van Krevelen diagram for different solid fuels (Note: The reference data of 
bituminous, sub-bituminous coal and wood are from IFRF solid fuel database). 


An ECO PHYSICS CLD 700 EL ht was used to measure the NOx, and 
the error of the measured NOx was less than 0.5% of the measuring 
range. Before each experiment, the test chamber was heated for 
10-12 h to achieve a uniform temperature profile inside the cham¬ 
ber. Two observation windows near the tip of the fuel lance al¬ 
lowed for the observation and recording of the flame behavior. 

2.3. Flame parameters 

The physical characteristics of the flame, which include geomet¬ 
rical and luminous profiles, temperature distributions and oscilla¬ 
tory attributes, provide instantaneous information on the quality 
of the fuel, and thus the quality of the combustion process [3]. A 
more accurate and less subjective method was proposed to de¬ 
scribe a flame based on the species concentrations; this type of 
flame was named as “chemical flame” [29]. The chemical flame 
of solid fuels and the flame parameters are schematized in Fig. 3. 
All of the flame parameters based on the chemical flame data, 
including the flame lift-off distance, flame length and maximum 
flame diameter, are distinguished from the luminous flame param¬ 
eters based on visual observations. In this work, the chemical flame 
was employed to describe the shape and position of a flame inside 
the furnace based on the following equation [29]: 




Fig. 2. Scheme and selected geometry of the test furnace. 
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Fuel/air injection port 



Fig. 3. Chemical flame scheme of solid fuels. 


Ro = m + m (1) 

m 0 + 2^i S ' m F,i 

where R 0 is the oxidation mixture ratio, m is the mass fraction, and 
the subscripts 0 and F refer to oxygen and the combustible species, 
respectively, s = n 0 M 0 /n F M F , in which n is the stoichiometric coef¬ 
ficient (number of moles) and M is the molecular weight. 

It was reported that the mechanism of pulverized coal flame 
propagation depends on the rate of appearance of combustible vol¬ 
atiles in front of the flame, and heterogeneous combustion is influ¬ 
enced by the rate of pyrolysis behind the flame[30]. Therefore, it is 
possible to employ the oxidation mixture ratio to describe the 
flame volume for pulverized fuel combustion. According to Yang 
et al. [29], the flame volume can be defined as 0 < R 0 ^ Rc and 
the critical oxidation ratio of R c = 0.99 approximately represents 
the outside border of the flame. 


2.4. CFD modeling 

Ansys Fluent 13.0 was used to simulate the combustion, fluid 
and particle flow, as well as the heat and mass transfer inside 
the furnace [31]. A cylindrical combustion chamber of 0.4 m inner 
diameter and 1.96 m height was modeled with a 3D geometry, as 
presented in Fig. 2. The mesh independence was verified by com¬ 
paring the reaction results. The mesh consisted of 73,588 hexahe- 
dral cells and was refined in the zone of fuel injection, where the 
ignition and most of the combustion reactions occurred. In this 
work, the standard k-s model with a standard wall function was 
used to calculate the effect of turbulent flow [32]. The discrete 


ordinates (DO) radiation model was chosen to simulate the radia¬ 
tion heat transfer. 

2.4.1. Combustion model 

The two-competing-rate model was employed to predict the 
biomass devolatilization rate and the volatile matter yield [33]. 
Even though the volatile matter contains many gaseous species, 
it was reasonable to represent the released volatile matter by a sin¬ 
gle virtual material, because the reaction rate of the volatile matter 
was limited by the turbulent mixing rate of the evolved gas and the 
oxidizer at the near-burner condition. 

In high temperature air combustion cases, the maximum fur¬ 
nace temperature is as high as 1400 °C, therefore, volatile combus¬ 
tion occurs at high Damkohler numbers (Da>l). The type of 
volatile matter, which mainly consists of gas phase hydrocarbons, 
depends on the fuel type and the pyrolysis conditions. Two-step 
global combustion is used for the homogeneous combustion of 
gaseous volatiles and involves two overall reactions that were cal¬ 
culated using Eddy-Dissipation model [31]. The char combustion 
was determined by the minimum kinetic rate and diffusion rate 
[34,35], and the diffusion coefficients for all cases were calculated 
using the following equation [36]: 

D 0 = 2.53 x lCr 7 r^ 75 /d p (2) 

2.4.2. Discrete phases model (DPM) 

Both biomass particles and coal particles were modeled using 
the DPM approach. The trajectory of the discrete particle was cal¬ 
culated by integrating the force balance on the particle as in the 
motion equation [31]: 

du P 18/i C D -Re g(Pp-p) m 

T Pr <?„ 24 + — ST“ <3 ' 

Based on Eq. (3), the drag force has an important effect on the 
trajectory of discrete particles. The drag force coefficient, C D , is a 
function of the particle shape factor <F> [37]. The shape factor is 
defined as the ratio of the surface area of a sphere with the same 
volume of the particle and the actual surface area. In this work, a 
non-spherical drag law is adopted to calculate the drag force 
coefficient of the biomass particles [37], and the shape factor of 
biomass particle was related to its volatile content, as presented 
in Table 2. For coal particles, it is widely acknowledged that the 
drag force can be accurately calculated using the spherical drag 
law [31,38]. 

Due to the significant fraction of volatile matter in biomass, the 
effect of swelling during the biomass devolatilization process on the 
drag force should also be considered. In this work, the swelling 


Table 2 

Kinetic parameters for torrefied biomass and coal particles. 



TB1 

TB2 

TB3 

Coal [39,45,46] 

Devolatilization kinetics 
Al (s- 1 ) 

(at high heating rates) 
6.40E+03 5.48E+03 

3.27E + 03 

3.7E+5 

A2 (s' 1 ) 

1.98E+05 

2.15E+07 

2.10E + 07 

1.46E+13 

El (kj/mol) 

57.0 

59.8 

60.2 

73.7 

E2 (kj/mol) 

145 

147 

159 

251 

VI 

0.67 

0.62 

0.499 

0.3 

V2 

0.74 

0.70 

0.50 

0.8 

Mean error 3 (%) 

4.5 

4.5 

4.5 

- 

Char oxidation kinetics 

A( s-V 1 

56 

56.3 

56.3 

0.91 

Ea (kj/mol )kJ/mol 

77 

78 

79 

103 

Shape factor 

0.53 

0.54 

0.58 

1 

Swelling coefficients 

1.0 

1.0 

1.8 

1.1 


a Calculated as the mean of absolute differences between model and experi¬ 
mental results of torrefied biomass from high heating rate tests in the IPFR: 
error = Y) n ]Xmod ~ Xexpl • 100%. 
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coefficient for each type of biomass was experimentally obtained, 
and was calculated as the ratio of the average diameter of the char 
particles (obtained at a high heating rate in the IPFR reactor) to the 
average diameter of the initial torrefied biomass particles, the mea¬ 
sured swelling coefficients of torrefied biomasses are presented in 
Table 2. The swelling coefficient of pulverized coal was selected as 
1.1 according to the experimental studies of Fletcher [39]. 

2.4.3. NOx model 

Thermal NOx refers to the NOx formed via the high temperature 
oxidation of the atmospheric nitrogen in the combustion air. It is 
widely acknowledged that the formation rate of thermal NO is 
determined by the extended Zeldovich mechanisms [40]. The pro¬ 
duction of NOx from nitrogen-containing fuels results from the 
conversion of fuel-bound nitrogen to NOx during combustion. Fuel 
nitrogen is usually assumed to react via the intermediates of HCN 
and NH 3 , which are finally either oxidized into NO or reduced into 
N 2 [41 ]. In this work, it was assumed that 90% of the volatile nitro¬ 
gen from coal will be converted to HCN, and the rest will form NH 3 
[42]. For the biomass particles, it was assumed that 10% of HCN 
and 90% of NH 3 were derived from the volatile nitrogen [43]. In 
addition, it was assumed that the nitrogen contained both in the 
coal char and the biomass char will be directly converted to NO 
due to the complex conversion mechanisms. Finally, prompt NOx 
was also considered in this work [44]. 

3. Results and discussions 

3.1. TG studies of torrefied biomasses 

The results of the devolatilization studies in the TG balance can 
be used to compare the reactivity of fuels. The weight loss curves of 
all of the biomass materials are compared in Fig. 4a. The volatiles 



100 200 300 400 500 600 

Temperature (°C) 

Fig. 4. Comparisons of torrefied masteries TG (a) and DTG (b). 


were released in the main weight loss step. As observed from the 
proximate analysis (see Table 1), the amount of released volatile 
matter reduced gradually from the raw PKS to TB1, TB2, and TB3. 
Furthermore, a slightly higher reactivity was observed in the raw 
PKS by considering the lower temperature at which the devolatil¬ 
ization begins. This was observed in the DTG curves of Fig. 4b, in 
which the onset temperature did not vary significantly for the torr¬ 
efied PKS. Generally, the earlier sub-peak of devolatilization was 
attributed to the decomposition of hemicellulose, which is the 
more reactive chemical component of the biomasses. Therefore, 
the delay of this sub-peak observed in the torrefied materials could 
be explained by a significant modification to the hemicelluloses 
structure during the torrefaction process. The main devolatiliza¬ 
tion peak was observed at approximately 360 °C, which corre¬ 
sponds to the decomposition of cellulose. A higher peak indicated 
a higher content of cellulose. The decomposition of lignin occurs in 
a wider temperature range, starting from 250 °C and continuing to 
the end of the process for all of the torrefied PKS materials. 

3.2. Kinetics of the torrefied biomasses 

In the two-competing-rate model, the first reaction (Al, El, and 
VI) is used to calculate the devolatilization rate at the lower tem¬ 
peratures, whereas the second reaction (A2, E2, and V2) plays a 
dominant role at higher temperatures. Generally, the volatile mat¬ 
ter yield is higher at high temperature than the standard volatile 
matter revealed by the proximate analysis. Once the pre-exponen¬ 
tial factor and activation energy are known, the model calculates 
the release rate and the yields of volatile matter. For the torrefied 
PKS particles, the parameters for the kinetic rate of devolatilization 
and the char oxidation were obtained from the experimental tests 
at high temperatures and high heating rates in the IPFR. The kinet¬ 
ics of the torrefied PKS are shown in Table 2. For the coal particles 
used in this study, the kinetic parameters of devolatilization were 
obtained from Li et al. [45], and the kinetics of the coal char oxida¬ 
tion were calculated using Hurt and Mitchell’s correlations [46]. 

3.3. Validation of the reference case 

The flames were photographed and video recorded using a CCD 
commercial photo camera (Sony DSC-V1) via two observation win¬ 
dows. In this paper, the observed flame is referred to as the visible 
flame. The flame tests were performed with the coal particles, and 
thus the visible coal flame will be used to validate the combustion 
model. 

To further understand the combustion performances when 
switching a pulverized coal-fired boiler to biomass operation, the 
combustion and flame characteristics of the torrefied biomass were 
studied numerically based on a laboratory-scale down-firing fur¬ 
nace under high temperature air conditions. In total 17 cases were 
simulated, and their detailed operating parameters are listed in 
Table 3. 

The temperature profile, visible flame, and chemical flame of 
pulverized coal were compared to validate the used model, as 
shown in Fig. 5, in which the temperature profile and the chemical 
flame were obtained based on the simulation results. The high 
temperature contour can approximately represent a luminous 
flame, and the high temperature zone is caused by the highly exo¬ 
thermic reaction that takes place. However, the flame cannot be 
described exactly based on the temperature profile, because the 
definition of the luminous flame boundary has always been chal¬ 
lenging. Compared to the temperature profile and the visible 
images of the flame, the chemical flame provides a more accurate 
description of the flame by considering of flame front, lift-off dis¬ 
tance, and maximum flame diameter. From Fig. 5, the flame lift- 
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Table 3 

The operation conditions of all cases. 


Case 

Fuel 

Fuel type 

Feed rates (kg/h) 

Particle size distribution (pm) 

Transport air 
Velocity (m/s) 

Temp. (°C) 

Secondary air 

Oxygen concentrations (vol%) 

Flow rates (kg/h) 

Temp. (°C) 

1 

Coal 

5.00 

R.-R. (avg. 61.7, max. 300) 

30 

120 

21/(Air) 

115 

1150 

2 

Coal 

5.00 

R.-R. (avg. 61.7, max. 300) 

20 

120 

21/(Air) 

115 

1150 

3 

Coal 

5.00 

R.-R. (avg. 61.7, max. 300) 

9 

120 

21/(Air) 

115 

1150 

4 

Coal 

5.00 

R.-R. (avg. 61.7, max. 300) 

20 

120 

17.9/(Air +Steam) 

115 

1150 

5 

Coal 

5.00 

R.-R. (avg. 61.7, max. 300) 

20 

120 

13.9/(Air +Steam) 

115 

1150 

6 

Coal 

5.00 

R.-R. (avg. 61.7, max. 300) 

20 

120 

11.8/(Air + Steam) 

115 

1150 

7 

TB1 

6.61 

R.-R. (avg. 61.7, max. 300) 

30 

90 

21/(Air) 

115 

1150 

8 

TB2 

6.73 

R.-R. (avg. 61.7, max. 300) 

30 

90 

21/(Air) 

115 

1150 

9 

TB3 

5.62 

R.-R. (avg. 61.7, max. 300) 

30 

90 

21/(Air) 

115 

1150 

10 

TB2 

6.73 

R.-R. (avg. 61.7, max. 300) 

20 

90 

21 /(Air) 

115 

1150 

11 

TB2 

6.73 

R.-R. (avg. 61.7, max. 300) 

9 

90 

21/(Air) 

115 

1150 

12 

TB2 

6.73 

Uniform-100 

20 

90 

21/(Air) 

115 

1150 

13 

TB2 

6.73 

Uniform-200 

20 

90 

21 /(Air) 

115 

1150 

14 

TB2 

6.73 

Uniform-300 

20 

90 

21/(Air) 

115 

1150 

15 

TB2 

6.73 

R.-R. (avg. 61.7, max. 300) 

20 

90 

17.9/(Air +Steam) 

115 

1150 

16 

TB2 

6.73 

R.-R. (avg. 61.7, max. 300) 

20 

90 

13.9/(Air + Steam) 

115 

1150 

17 

TB2 

6.73 

R.-R. (avg. 61.7, max. 300) 

20 

90 

11.8/(Air + Steam) 

115 

1150 


Note: R.-R. represents Rosin-Rammler size distribution with spread parameter of 4.21. 


Temperature [K) 

1 2145 
2034 
1923 
j- 1812 
\ 1702 
r 1591 
r 1480 
r 1370 
r 1259 
r 1148 
r 1037 
r 927 

j- 816 

r 705 



Temperature profile 


Chemical flame 


Fig. 5. Comparisons among chemical flame, temperature profile and visible flame. 



off distance and flame profile of coal powders were similar when 
the chemical flame was compared with the visible flame. 

Furthermore, the lift-off distances of the chemical flame and the 
visible flame were compared at various transport velocities, as 
shown in Fig. 6. It can be observed that the lift-off distances of 
the chemical flame were in good agreement with that of the visible 
flame, and the flame lift-off distances progressively increased as 
the transport velocity increased. A larger difference was observed 
when the pulverized coal was injected by the transport air at a 
velocity of 30 m/s, which might be explained by the increase in 
the measured error, because a high flame propagation speed en¬ 
hanced the fluctuation of the visible flame. Based on the above dis¬ 
cussions, it could be concluded that chemical flame based on 
numerical results was reasonable for further studies of biomass 
flames. 


3.4. Effect of the torrefaction degree on the flame 

Generally, hemicellulose is decomposed during the biomass tor- 
refaction process, and thus a high torrefaction degree usually refers 
to a small amount of volatile in remaining biomass. Four different 
types of solid fuel were studied: coal and three torrefied biomasses 
(TB1, TB2, and TB3) with different torrefaction degrees. The four 
cases (cases 1, 7, 8, and 9) were simulated with the same heat fur¬ 
nace capacity, transport air velocity, and particle size distribution. 
Figure 7 shows the effects of the fuel type on the flame position 
along the furnace height and the flame shape. It should be noted 
that the 0 mm height of the down-firing furnace is the fuel/air 
injection ports. Figure 7 indicates that the coal flame was closest 
to the fuel inlets, whereas all three torrefied flames moved to posi¬ 
tions far from the fuel inlet, and the distance that the biomass 
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Transport air velocity (m/s) 


Fig. 6. Comparison of chemical flame lift-off distance and visible flame lift-off 
distance. 


can effectively surround the coal particle against the oxygen attack. 
According to Krazinski et al. [47], the volatiles are released from 
the surface of the particles and subsequently combusted with a dif¬ 
fusion flame, and the volatiles are released rapidly in the front of 
the flame; in contrast, the devolatilization rate is no longer signif¬ 
icant in the tail of flame, which may also affect the flame lift-off 
distance. Meanwhile, it should be noticed that the predicted peak 
flame temperature of the coal powders was higher than that of 
the biomass powders, and the peak temperatures of the torrefied 
biomasses decreased as the volatile contents increased. The pre¬ 
dicted peak temperatures of coal, TB3, TB2, and TB1 under the 
studied conditions is 2186 °C, 2007 °C, 1963 °C, and 1843 °C, 
respectively. The temperature of flame formed around the particles 
depends on the heating value of gas releasing from particles. 
Therefore, the delay of the high temperature flame zone, the 
changes in the flame volumes with volatile content of the fuel, 
and the relatively lower peak flame temperature must be taken 
into considered when combusting biomass in an existing pulver¬ 
ized coal boiler. 


flame moved was approximately a length of the coal flame, which 
is confirmed by the experimental results reported by Lu et al. [3]. 
The results might be caused by the differences in the drag forces 
imposed on the biomass particles and coal particles. Due to the 
lower density of biomass, an enhanced drag force on the biomass 
particles could cause a late release of volatile matter and thus delay 
the ignition of the fuel-air mixture. This conclusion further ex¬ 
plains that the major combustion zone of the biomass that moved 
towards the furnace exit in a pulverized co-firing boiler, as re¬ 
ported by Li et al. [10]. 

Moreover, it was found that the torrefied biomass flame had a 
wider diameter than coal flame, a similar experimental result 
was reported by Ballester et al. [12]. It is interesting to note that, 
Fig. 7 shows the biomass flame volume increased as the volatile 
content increased for all of the studied torrefied biomasses. A clear 
reason is that a higher volatile content improved the diffusion of 
the combustible gases, causing a larger flame volume and a large 
flame diameter; in addition, volatile and its preliminary combus¬ 
tion products enhances the gas phase diffusion from near-flame 
zone to a wider zone. At high rates of devolatilization, the volatiles 


3.5. Effect of the oxygen concentration of the oxidizer on flame 

Volumetric combustion is characterized as an extension of the 
HTAC technology for solid fuels with typical advantages, including 
the dilution of the oxygen by an enhanced flue gas internal recircu¬ 
lation, uniformity of both the temperature and species concentra¬ 
tions, and a lower peak flame temperature. In this work, to 
understand volumetric combustion properties of torrefied biomass 
particles, the flame characteristics of a torrefied biomass (TB2) 
were studied numerically with varying the oxygen concentrations 
(cases 10, 15, 16, and 17). The oxygen concentration was adjusted 
by using a mixture of steam and air. The effect of oxygen concen¬ 
tration in the oxidizer on the flame position along the furnace 
height and the flame shape is shown in Fig. 8. The flame volume 
and flame length increased as the oxygen concentration in oxidizer 
was reduced, whereas the flame lift-off distance and maximum 
flame diameter changed only slightly for all of the studied cases. 
This is because the oxygen concentrations were adjusted by replac¬ 
ing part of the air with high temperature steam, and the combus¬ 
tible species released after devolatilization diffused in a wider zone 


0.1 v!v! - Fuel/air inlets at 0 m - 
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Fig. 7. Flame shapes and flame positions along furnace height with various fuel F ig. 8. Flame shapes of TB2 and flame positions along the furnace height with 
types- various oxygen concentrations. 
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of the furnace before being completely combusted. Therefore, the 
combustion volume was larger for the cases with diluted oxygen 
concentrations than that with normal air combustion. In addition, 
according to the simulation results, the peak temperatures of bio¬ 
mass decreased as expected with deceasing oxygen concentrations, 
which is a common feature of volumetric combustion with low 
oxygen concentrations. The predicted peak temperature of TB2 at 
oxygen fractions of 21%, 17.9%, 13.9% and 11.8% is 1913 °C, 
1649 °C, 1387 °C, and 1220 °C, respectively. 

The benefits of low peak temperature and low oxygen concen¬ 
tration have been demonstrated in various industrial combustion 
applications, because it is an effective option to reduce NOx and 
improve boiler efficiency. To understand the relationship between 
the flame length or the flame volume and the NOx emissions, the 
performances of the NOx emissions under different oxygen con¬ 
centration conditions have also been studied. Figure 9 indicates 
that the oxygen concentration strongly influences the NO emission 
and the flame length, and a lower oxygen concentration always 
corresponds to a longer flame length and a lower NO emission in 
studied biomass and coal cases. According to the predicted results, 
the lowest amounts of NO produced are 248 mg/N m 3 for biomass 
and 515 mg/N m 3 for coal when the oxidizer contained 11.8% oxy¬ 
gen. Accordingly, all of the measured NOx emissions from coal 
showed a similar trend with the predicted values, which further 
confirmed the reasonableness of the numerical models and meth¬ 
od adopted in this study. Although Fig. 9 shows that the maximum 
amount of NOx was higher than 1400 mg/N m 3 in case of normal 
air combustion, was mainly due to the injection of the fuel and 
the oxidizer without applying any low-NOx technology, such as 
air-staging and fuel-staging. 

According to the results, it could be concluded that a longer 
flame length or a bigger flame volume is beneficial for NOx reduc¬ 
tion. There might be several reasons to explain these results. The 
first important reason is that the whole combustion zone was 
operated under sub-stoichiometric conditions with an excess air 
ratio less than one, which not only reduced the maximum flame 
temperature, but also prevented the formation of nitrogen oxides 
from fuel-bound nitrogen. Furthermore, the lower 0 2 partial pres¬ 
sure and lower maximum flame temperature under lower oxygen 
combustion conditions inhibited the formation of thermal NO. 
Third, because the triatomic species H 2 0 is not transparent to radi¬ 
ation, it substantially increased the absorptivity and emissivity of 
the flue gas, and as a result, an increased H 2 0 concentration may 
have affected the heat transfer and flame temperature inside the 
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Fig. 9. Effects of oxygen concentrations on NO emissions and flame lengths. 


furnace, which could in turn affect the formation of NO [48]. The 
results further confirmed that volumetric combustion with a flame 
as large as furnace volume could be featured as a promising low- 
NOx combustion technology for biomass combustion. 


3.6. Effect of the transport air velocity on flame 

The transport air velocity is important in the combustion pro¬ 
cesses of coal and biomass powders, and can partly control the glo¬ 
bal combustion rates. In this work, a torrefied biomass (TB2) was 
studied. Three cases with different transport air velocities of 9 m/ 
s, 20 m/s, and 30 m/s (cases 9, 10, and 11) were simulated at the 
same fuel feed rate and same particle size distribution. The trans¬ 
port air velocities were varied by changing the amounts of the air 
flow rates. The changes of flame position along the furnace height 
and flame shape with varying transport air velocities are shown in 
Fig. 10. Obviously, the transport air velocity had significant effects 
on both the flame volume and the flame length, and a longer flame 
was observed at a lower transport air velocity. Under the high tem¬ 
perature combustion conditions, the operating temperature of the 
furnace is high (above 1400 °C), so that volatile combustion occurs 
at high Damkohler numbers (Da > 1); therefore, the combustion 
rates of the volatiles are mainly controlled by the diffusion of oxy¬ 
gen. At lower transport air velocities, the diffusion rate of oxygen is 
slower than that under a higher transport air velocity. However, 
the flame front decreased as the transport air velocities increased, 
and the flame lift-off distance increased simultaneously. Further¬ 
more, the transport air velocity had negligible effects on the max¬ 
imum flame diameter because the maximum flame diameters 
were almost the same for the three cases, as shown in Fig. 10, 
which might further confirm that the maximum flame diameter 
is mainly affected by the volatile contents of the fuels. In addition, 
according to the numerical simulation, the predicted peak temper¬ 
atures of TB2 at transport velocities of 9 m/s, 20 m/s, and 30 m/s is 
1855 °C, 1913 °C, and 1963 °C, respectively. A possible explanation 
is that the high transport air velocity enhanced the diffusion of the 
oxidizer onto the char surface and resulted in a high char oxidation 
rate inside the flame. According to the above discussion, a rela¬ 
tively lower fuel injection velocity is suggested when combusting 
biomass in an existing pulverized coal furnace, meanwhile, it indi- 
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Fig. 10. Flame shapes for TB2 and flame positions along furnace height with various 
primary air velocities. 
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cates that the discrete concentration in pneumatic transport flow 
would be increased to maintain the same heat capacity. However, 
the residence time, combustion rate, furnace temperature and eco¬ 
nomical factors should be considered altogether to determine a 
reasonable transport air velocity; for example, a lower flow rate 
usually means a lower combustion rate, whereas a high flow rate 
normally accompanies a high erosion rate. 

3.7. Effect of particle size on flame 

Particle size is an important parameter during the pulverized 
fuel combustion process, and thus the effects of particle size on 
combustion performances have been studied widely. Typically, it 
is generally acknowledged that heterogeneous reactions are con¬ 
trolled by diffusion for particles of diameter larger than 0.1 mm 
when surrounding temperatures are higher than 1000 K; while 
the heterogeneous reactions of even 0.01 mm sized particles are 
controlled by diffusion when the surrounding temperature exceeds 
2000 K [47]. 

To investigate the influence of particle size on high temperature 
flames, three cases (cases 12, 13, and 14) have been studied with 
three different uniform particle sizes (0.1 mm, 0.2 mm, and 
0.3 mm, respectively). Larger biomass particles were not investi¬ 
gated in this work because of the limitation of the fuel lance diam¬ 
eter. The effect of particle size on the flame position along the 
furnace height and the flame shape is shown in Fig. 11. It can be 
observed that the flame positions gradually moved toward the fur¬ 
nace exit as particle size increased. This result could be explained 
by the longer residence time needed for the complete pyrolysis 
of a larger particle, which would thus delay the burnout of the 
char. For each type of solid fuel, the efflux rate of the volatile 
was found to vary with the particle sizes. For example, Howard 
and Essenhigh reported that coal particles larger than 65 pm do 
not react heterogeneously during the rapid devolatilization process 
[49]. Furthermore, according to the predicted results, the peak 
temperature for TB2 consisting of uniform particle sizes of 
0.1 mm, 0.2 mm, and 0.3 mm is 1984 °C, 1881 °C, and 1582 °C, 
respectively, which indicates that a small particle is quickly com¬ 
busted completely. In addition, the results showed that the maxi¬ 
mum flame diameter was almost the same for the studied 
particle sizes; therefore, it could be concluded that particle size 
had negligible effects on the maximum flame diameter. 



Fig. 11. Flame shapes for TB2 and flame positions along the furnace height with 
various particle sizes. 


4. Conclusions 

In this work, the combustion and flame characteristics of torr¬ 
efied biomass have been studied numerically in a laboratory-scale 
down-firing furnace under higher temperature air conditions. 
Three torrefied biomasses obtained from same raw PKS material 
with torrefaction degrees of 10%, 20% and 30% of volatile content 
released on a dry basis were studied. Finally, the effects of the tor- 
refaction degree, oxygen concentration, transport air velocity, and 
particle size on the flame positions and flame shapes, and peak 
flame temperatures were discussed. 

The results mainly showed that the enhanced drag force on the 
biomass particle caused a late release of volatile matter and re¬ 
sulted in a delay in the ignition of the fuel-air mixture, and the 
maximum flame diameter was mainly affected by the volatile con¬ 
tent of materials. Furthermore, a lower oxygen concentration al¬ 
ways corresponded to a larger flame volume, a lower peak 
temperature, and a lower NO emission. A longer flame was ob¬ 
served when the transport air velocity was lower, and the flame 
positions gradually moved toward the furnace exit with increasing 
of particle sizes. The results could be used to design a biomass 
combustion chamber and/or provide an operational reference for 
the co-firing performances in a coal-fired boiler. However, a series 
of relative experiments also should be carried out in the future to 
obtain further details of the biomass flames. 
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